Spatial stochastic resonance in ID Ising systems 
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The ID Ising model is analytically studied in a spatially periodic and oscillatory external magnetic 
field using the transfer-matrix method. For low enough magnetic field intensities the correlation 
between the external magnetic field and the response in magnetization presents a maximum for a 
given temperature. The phenomenon can be interpreted as a resonance phenomenon induced by 
the stochastic heatbath. This novel "spatial stochastic resonance" has a different origin than the 
classical stochastic resonance phenomenon. 
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I. INTRODUCTION 

Many recent papers PH|] revealed une quiv ocally the 
phenomenon of stochastic resonance (SR) [1(J] in the ki- 
netic Ising model driven by a temporary oscillating mag- 
netic field. SR was anticipated by considering the Ising 
model as a system of coupled two-state oscillators in a 
stochastic force-field which is taken to be thermal fluc- 
tuations. In this sense the system has all the ingredients 
necessary to observe the classical phenomenon of stochas- 
tic resonance. 

In the present paper we intend to study the ID fer- 
romagnetic Ising model in a spatially periodic and os- 
cillatory B(i) external magnetic field (B(i + A) = B(i), 
< B(i) >i— 0, 2A the spatial period of the field, and i 

the lattice points i = 1,2,3, 2p\, p — 1,2, ). We 

impose periodic boundary conditions, thus S(2p\ + 1) = 
5(1). The Hamiltonian of the system writes as 



H 



2p\ 

-jJ2s(i)S(i 



1) 



2pA 

»Y2B(i)s(i 



(i) 



(p the magnetic moment of the S(i) = ±1 Ising spins). 
The magnetic field is taken stationary in time. 

For the proposed problem is easy to realize that at 
T = thermodynamic temperature and not to high mag- 
netic field intensities, due to the oscillatory nature of the 
external magnetic field the a =< B(i)S(i) > (the brack- 
ets denotes both a spatial and ensemble average) correla- 
tion is greatly reduced, i.e. the infinite correlation length 
(£(0) = °°) between the spins competes with the finite 
2A period of B(i). For T — oo is also obvious that due 



to the leading stochastic contribution a = 0. We ex- 
pect that for a given finite temperature the £ correlation 
length will be of the order of the A period of B(i) and 
thus the a correlation will reach a maximal value. This 
resonance like phenomenon is induced by the stochastic 
force field (temperature) for the energetically frustrated 
system at T = 0. Therefore, in order to distinguish 
this phenomenon from the classical SR phenomenon we 
named it spatial stochastic resonance . 



II. THE METHOD 

To give an exact solution for the proposed problem 
we choose the most simplest possible B{i) configuration 
with the above imposed properties. We choose B(i) = B 
for i = 2n\ + 1 (n = 0, 1,2,. ..p - 1), B(i) = -B for 
i = (2n + 1)A + 1, and B(i)=0 for all other lattice points. 
We are interested in the < 5(1) > average magnetization 
at the i — 1 position from where the a correlation is easily 
determined. From the chosen magnetization profile we 
get: 

a =< B(i)S(i) >= P B(< 5(1) > - < 5(A + 1) >) (2) 

<5(A + 1)>, 



From symmetry arguments < 5(1) > 
and we can write: 

a p = 2 P B < 5(1) > 



(3) 



In order to determine < 5(1) > we calculate i.) the 
Z2p\ partition function of the system ii.) the Z^. ' x par- 
tition function for 5(1) = 1 imposed condition, iii.) and 



the Z 2 x partition function for the 5(1) = — 1 imposed 
condition. We get the desired < 5(1) > value, as: 



< 5(1) >= 



^2p\ Z 2p \ 
Z 2p \ 



(4) 



During our calculations we use several matrices whose 
explicit forms are given in the Appendix. 



III. THE L = 2A LENGTH CHAIN (P = 1) 

With the notations j = J/kT and h = ^iB/kT (T the 
temperature of the system, k the Boltzmann constant) 
the partition function Z 2 \ writes as 



S(1)S(2) S(2A) 



S(i)S{i+l)]+hS(l)-hS(X+l) 



(5) 



(the sums are for S(i) = ±1). 
We use the transfer-matrix method to calculate Z2\ 

Z 2 \ = Es(i) Es(2) •■■ Es(2A) < s ( l ) I 7 + I s ( 2 ) > 

< 5(2) | I | 5(3) > ... < 5(A) | I | 5(A + 1) > 

< 5(A + 1) | J_ | 5(A + 2) >< 5(A + 2) | I | 5(A + 3) > 

... < 5(2A - 1 | I | 5(2A) >< 5(2A) | I | 5(1) > 



with the I and I± matrices given in the Appendix (32) 
It is now immediate to realize that, 



^2A 



Tr(I + lfl-Io) 



(G) 



(where Tr(A) denotes the trace of the A matrix). Intro- 
ducing the M (|33|) diagonal matrix we get 



(7) 

(8) 

(9) 

(10) 



1+ = MI 

I_ = M- l I Q 

Z 2 x = Tr(W) 

W = MI X M^I X 



Choosing a representation where I becomes diagonal we 
get: 



I' = UIoU- 1 

M' = UMXJ- 1 

W' = M'l'^M'- 1 !'^ 

Z 2 x = Tr(W) 



(11) 
(12) 

(13) 
(14) 



(the U, I' Q and M' matrices are also given in the Ap- 
pendix). After some elementary algebra one will find: 



Z 2 \ = 2 2X {cosh 2 (h)[cosh 2X (j) + sinh 2X (j)} 
— 2sinh 2 (h)sinh x (j)cosh x (j)}} 



(15) 



For Z 2X and Z 2X 



z ± = y^y ... y^ e jE-!i s w 5 ( i+i )]+ ,iS ( i )- ,iS ( A + i ) 

5(2)3(3) S(2A) 

(16) 

(for the ± cases we have 5(1) = ±1 respectively) we 
perform a similar calculation: 

^2A = Es(2) Es(3) — Ss(2A) < S(2\) | S± | 5(2) > 

< 5(2) I I I 5(3) > ... < 5(A) | I | 5(A + 1) > 

< 5(A + 1) | I_ | S(X + 2) >< 5(A + 3) | I | 5(A + 4) > 

... <5(2A-1) \I | 5(2A) > 



J 2A 



Tr(S ± I x - l I-i;- 2 ) 



(17) 



The S± matrices are also given in the Appendix (J3C 
Using (pj) we can write: 



Zf x = Tr{P ± I x M-'I x ) 



p± = i- x s ± i- x 



(18) 
(19) 



Again, we calculate the trace in the representation where 
I is diagonal and we get: 



Z 2\ = 2 2X - l e ±h {cosh{h)[cosh 2X (j) + sinh 2X (j)} =p 
2sinh(h)sinh x (j)cosh x (j)} 



(20) 



With the obtained Z 2 \ and Z 2X values < 5(1) > is eas- 
ily determined (Q) and for the a v @ correlation we get 
finally: 



o-i = 2B 



tanh(2h) 



*- ' cosh(2h) V 1 



l+tanh x (j) 



(2/i) \l-tanh x (j) 



(21) 



IV. THE L = 2PX (P > 1) LENGTH CHAIN 

In order to be able to make the calculations easily in 
the p > 1 case we first write Z 2X in a more convenient 
form. 



Z= x = Tr{P'±r x M'- x r x ) = Tr{R'±W) (22) 

(23) 



R'± = P^M'- 1 



Is easy to realize that: 

Z 2X = Z+ x + Z 2X = Tr(W) 



(24) 
(25) 



ZJx - Z 2X = Tr(R'W') 

For the p > 1 case is immediate that: 

^2 P a = Tr(W'P) (26) 

Writing up the effective forms of Z 2 x like in the p = 1 
case, one can also show that: 



Ztx = Tr(R' ± W'n 



(27) 



In the representation where W is diagonal it is easy now 
to calculate < 5(1) > and a p . Denoting by xi an d X2 
(Xi > X2) the eigenvalues of W, after some simple alge- 
bra we find: 



L P D Tr{W'P) 



G\ 



Tr(W') X\~xl 



VA x\+xl 



(28) 

A = Tr(W') 2 - 4Det(W) (29) 

Xi, 2 = T ^ w '^ (30) 

In the limit p — > 00 we get the simple formula 

Tr{W) 



<7i- 



V(A) 



(31) 



which is easily computable from the ( pTj ) VI" matrix given 
in the Appendix. 



V. DISCUSSION 

Equation ([H]) and (|l|) gives us the cr =< B(i)S(i) > 
correlations for the p = 1 and p = 00 periodic chains. 
The c(T) correlation is proportional with the < 5(1) > 
(T) curves (|). On Fig. 1 we plotted < 5(1) > (T) for 
three different applied magnetic field intensities. 

As expected, for fj,B/J < 2 (when the interaction with 
the external field is weaker than the interaction with the 
neighboring spins) a clear resonance like behavior is ob- 
tained. Both for p = 1 and p = 00 < 5(1) > (T) exhibits 
a clear maximum at a T r 7^ resonance temperature. 
It is also observable that the < 5(1) > (T) curves for 
p = 1 and p = 00 are very close, thus the p = 1 result is 
qualitatively well describing the p > 1 cases as well. The 
T r resonance temperature depends both on the B inten- 
sity of the applied magnetic field and the characteristic 
A distance of the spatial oscillations of B{i). On Fig. 2 
we illustrate the T r (B) dependence, and on Fig. 3 the 
T r (X) trend. From Fig. 2 we learn that in the B — > 
limit the T r values are converging to a constant (which is 
dependent on p), and in the fj,B/ J > 2 limit T r = 0, thus 
no resonance behavior is obtained. The T r (X) variations 
(Fig. 3) are also the one expected from our phenomeno- 
logical considerations. In the limit A — > 00 we get T r — > 0, 
and T r is monotonically decreasing with increasing A val- 
ues. It is interesting to note that for the minimal possible 
A value (A = 1) the resonance like behavior is still present 
but the p = 1 and p — 00 curves are much more distant 
comparable with large A values case. 



VI. CONCLUSIONS 

In the present paper we studied the response of ID 
Ising chain to a spatially periodic and oscillatory mag- 
netic field. Considering the most simple magnetic field 



profile we exactly solved the problem by using the 
transfer-matrix method. We found that the a =< 
B(i)S(i) > correlation between the applied magnetic 
field and the local magnetization exhibits a maximum for 
a given T r resonance temperature (Fig. 1). The T r res- 
onance temperature depends monotonically on the spa- 
tial oscillation length of the magnetic field (Fig. 3). The 
value of T r depends also on the B intensity of the mag- 
netic field, and becomes independent of B in the small B 
values limit (Fig. 2). The length of the chain L = 2pX, 
characterized by the p integer value has small influence on 
the observed resonance-like behavior for large A values. 
The obtained resonance-like behavior is induced by the 
stochastic heat-bath and has a different origin from the 
classical stochastic resonance phenomenon. We named it 
spatial stochastic resonance. 
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e -JTh 
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\e h 


n" 




r 1 


1 1 




M = 





e- h _ 


> U =T2 


1 


-1 


) 



I' =2 



M' = 

M'- 1 = 

cosh(j) 




cosh(h) sinh(h) 
sinh(h) cosh(h) 

cosh(h) —sinh(h) 
—sinh(h) cosh(h) 



sinh(j) 



;5± 



= e ± h 



,±2j 



1 



W' = 2 



2A 



11 



12 



»'■; 



(32) 
(33) 
(34) 
(35) 
(36) 
(37) 



w' n = cosh 2 (h)cosh (j) — sinhr (h)cosh (J)sinh A (J) 
w[ 2 = sinh(h)cosh(h)sinh x (j)[sinh x (j) — cosh x (j) 
w'21 = sinh(h)cosh(h)cosh x (j)[cosh x (j) — sinh x (j)] 

w'22 = cosh 2 (h)sinh 2X (j) — sinh 2 (h)cosh x (j)sinh x (j) 

1 ±1 



R' = 



n = 


2 


'0 1" 

1 


; R 



±1 1 



1 



1 ±1 
±1 1 



(38) 
(39) 
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Figures Captions: 

Fig. 1 Characteristic shape of < S(l) > (T) for three 
different magnetic field intensities (iiB/k = 0.1; 1.0; 2.1). 
The continuous lines are for p = 1, the dashed ones for 
p = oo, A — 20 for all curves. 

Fig. 2 The T r resonance temperature as a function of 
the applied magnetic field intensity. We draw the results 
for two different A values. The continuous lines are for 
p = 1, the dashed ones for p = oo. 

Fig. 3 The T r resonance temperature as a function 
of the A length. We draw the results for two different 
applied magnetic field intensities. The continuous lines 
are for p — 1, the dashed ones for p = oo. 
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FIG. 1. Characteristic shape of < S(l) > (T) for three different magnetic field intensities (fiB/k = O.f ; 1.0; 2.1). The 
continuous lines are for p = 1, the dashed ones for p = oo, A = 20 for all curves. 
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FIG. 2. The Tr resonance temperature as a function of the applied magnetic field intensity. We draw the results for two 
different A values. The continuous lines are for p — 1, the dashed ones for p = oo. 




FIG. 3. The T r resonance temperature as a function of the A length. We present results for two different applied magnetic 
field intensities. The continuous lines are for p = 1, the dashed ones for p = oo. 



